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Abstract 

The COBE satellite's maps of the cosmic microwave background (CMB) anisotropy do not have the resolu- 
tion to discriminate between theories of cosmic structure formation based on inflation, and those based on field 
ordering following a symmetry breaking phase transition. For this purpose it is critical to resolve the CMB 
anisotropy on degree scales. In this paper we report on detailed calculations of the degree scale anisotropies 
predicted in the global string, monopole, texture and nontopological texture scenarios of structure formation, 
emphasising their distinct character from those predicted by inflation, and commenting on the prospects of 
their detection in the near future. 
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I. Introduction 

The COBE satellite's detection of large angular scale anisotropy in the cos- 
mic microwave background [1] has opened a new window on the early universe. 
But the picture COBE itself provides is limited by a large angular smoothing 
scale and a low signal to noise ratio. It falls short of discriminating between a 
large number of theories of structure formation, based on very different physics. 
For example, theories based on inflation and on cosmic fleld ordering are equally 
consistent with COBE's findings. Both sets of theories predict an approxi- 
mately scale invariant spectrum of multipole moments, and on COBE scales a 
very Gaussian anisotropy pattern. In particular, the distinctive nonGaussian 
signatures expected in the field ordering theories (lines for strings [2] , spots for 
textures [3]) would be smeared out by the COBE beam and instrument noise, 
and very hard to distinguish from the Gaussian noise pattern expected from 
inflation [4] [5]. For this purpose, resolution on an angular scale smaller than 
that subtended by the horizon at last scattering appears to be essential. 

None of the preexisting theories can claim the COBE findings as an un- 
mitigated success. Roughly speaking, the COBE result was a factor of two 
higher than expected in the simplest inflationary theory, and a factor of two 
lower than expected in the simplest field ordering theories. These theories have 
O = 1, A = 0, with one free parameter determining the level of fluctuations - if 
this is fixed from COBE, the 'standard' inflation plus cold dark matter (CDM) 
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theory predicts galaxy cluster velocity dispersions much higher than observec 
[6]. It is fairly easy to fix the situation by adding new adjustable parameters 
In the inflationary case, considering 'mixed' dark matter [7], inflationary po 
tentials which yield an important additional gravity wave component [8] , [9] , o 
adding a cosmological constant [10], all improve the situation. 

The simplest field ordering theories are in principle as predictive as th( 
simplest inflationary scenario, since the level of mass fluctuations and of mi 
crowave anisotropies depend on a single free parameter, the symmetry breakinj 
scale An appealing aspect of these scenarios is that one naturally produce; 
the correct amplitude of density fluctuations with of order the grand unifica 
tion scale, around 10^® GeV. But unlike the simplest inflationary theory, thesi 
theories rely on nonlinear effects in an important way, and are harder to calcu 
late with. The most detailed calculations so far suggest that if the paramete 
(po is normalised to fit the COBE 10" variance, the level of mass fluctuation; 
inferred on intermediate scales (~ 20/i~^ Mpc) is too small to explain the ob 
served galaxy streaming motions [5]. It has been suggested that this problen 
would be ameliorated in an open universe [11], by suppressing the anisotropic 
produced by the field ordering at late times and thus allowing a larger valui 
for A positive cosmological constant could be expected to have the sam( 
effect. 

The pattern of CMB anisotropy on degree scales offers an independen 
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and in principle far cleaner test, being unaffected by the complex processes 
involved in galaxy formation. The simplest measure of anisotropy, and the 
only relevant quantity in Gaussian theories, is the power spectrum of multipole 
moments. If the standard inflation plus CDM theory is correct, there should 
be a peak in the spectrum at I ^ 200 produced at the epoch of recombination, 
Z ~ 1000, due in part to the Doppler effect of photons scattering off moving 
electrons. In the theories considered here which involve topological defects 
(-/V = 2, 3, 4), reionisation is likely, and last scattering would be shifted forward 
to Z ~ 100. The angular scale subtended by the horizon at last scattering is 
(for n = 1, A = 0) = (1 + Zis)~i, of order 2° without reionisation and 6° 
with reionisation. As will be seen, the spectrum predicted by the field ordering 
theories with the assumption of reionisation, is quite distinct, with much less 
power in higher multipoles. 

Even more distinctive is the nonGaussianity in the field ordering theories. 
The assumption that the anisotropy pattern is in the form of random Gaussian 
noise is frequently made. However, if the microwave sky is nonGaussian, there 
will be additional to be gleaned from anisotropy maps. In the nonGaussian 
theories explored here for example, details of the pattern would tell us about 
the pattern of symmetry breaking in particle physics at very high energies, 
information presently available from no other source (see e.g. [12]). 

The data gathered so far from experiments at the South pole present an 
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incomplete but nonetheless intriguing picture. If the the upper limit {5T/T) < 
1.4 X 10~^ derived from one of the four channels of the Gaier et. al. experimen 
[13] is correct, it is enough to exclude the standard inflationary spectrum o 
fluctuations (normalised to COBE) with either hot or cold dark matter at thi 
95 % confidence level [14]. And the thirteen point scan from the Schuste 
et. al. experiment [15] contains an apparently nonGaussian 'bump'. Similarly 
data from the MAX experiment seem inconsistent with with Gaussian statistics 
Under that assumption discrepant upper and lower 95 % confidence bounds wen 
obtained on the level of CMB anisotropy at from two different parts of thi 
sky [16]. It is clearly of interest to determine the predictions of the nonGaussiai 
field ordering theories on these angular scales, while the observational situatioi 
is still being resolved. 

One of the sharpest distinctions between the Gaussian and nonGaussiai 
scenarios mentioned above is the likelihood that the universe was reionised afte: 
recombination. If around 0.1 — 1.0 per cent of the mass in the universe underwen 
nonlinear gravitational clustering and star formation at redshifts greater than o 
the order 100, this would likely release sufficient ionising radiation to reionise thi 
intergalactic medium, causing the photons to rescatter. In this case the epoch o 
last scattering would be shifted forward to a redshift Zig ~ 100(.05//f2B/i5o)^ 
where / is the fractional ionisation (assumed not too different from unity). Ii 
field ordering theories with topological defects like strings or textures, large loca 
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density perturbations are produced at high redshift, and complete reionisation 
(/ = 1) is possible [5]. But in Gaussian scale invariant theories of structure 
formation, reionisation is very unlikely because high amplitude fluctuations are 
exponentially suppressed - using the Press-Schechter [17] or number of peaks 
formulae [18] one deduces that for the standard CDM theory normalised to 
COBE, no more than of order 10"^*^ of the baryons in the universe could have 
undergone collapse onto dark matter lumps, and thus be bound into stars, by 
a redshift of 100. This is far too little to reionise the universe. 

In this paper we shall make the assumption that the universe was reionised 
at redshifts greater than Zis given above, and study the consequences for the 
pattern of microwave anisotropy. We shall give results for the simplest sym- 
metry breaking theories: where an assumed 0{N^ global symmetry is spon- 
taneously broken by an N component scalar field taking a GUT scale vacuum 
expectation value. For N = 2, this is the theory of global cosmic strings, N = 3, 
global monopoles, N = 4:, global texture and N = 6, a,n example of 'nontopo- 
logical' texture. Of these the 'generic' cases may be argued to be global strings, 
which result from the complete spontaneous breaking of any abelian continu- 
ous global symmetry, and texture, which results from the complete breaking of 
any nonabelian continuous global symmetry. It is important to keep in mind, 
however, that 'realistic' unified theories may include more complex symmetry 
breaking schemes, in which results different by a factor of two might well be 
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obtained. For example the family symmetry scheme studied in [12] correspond, 
to two N = 6 theories, coupled in such a way as to produce topological texture 

In previous work, Bouchet, Bennett and Stebbins [19] constructed ai 
anisotropy map for gauged cosmic strings. They used a flat spacetime ap 
proximation to compute the anisotropy and ignored electron scattering, whicl 
in a reionised universe acts to smear out the distinctive linelike discontinu 
ities produced by cosmic strings [20], for redshifts earlier than last scattering 
Reionisation seems likely in the string scenario: in this case the strings wouk 
only start to produce clear discontinuities well after last scattering, on scale; 
in excess of 9is ~ 6". And of course these discontinuities would be very heavil; 
masked by the smaller angular scale fluctuations produced at earlier epochs. 

In our calculation for N = 2,we do not aim at infinite resolution, but rathe 
to produce a map smoothed on a scale comparable to the beam size adoptee 
for degree scale measurements. We include the effects of electron scattering 
ignored in the earlier calculations, and we believe that this provides a mon 
realistic picture of the anisotropy pattern one could expect from a degree scaL 
experiment. 

Techniques 

We study the evolution of the symmetry breaking scalar fields using, fo 
N > 2, the finite difference scheme of [5] for the nonlinear sigma model. W 
are able, with these techniques, to evolve the scalar fields in boxes in excess o 
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200'^ grid points, and for A'' > 2 to convincingly demonstrate that the results 
are independent of box size. For N = 2, we have developed a more accurate 
difference scheme to evolve a single 'angular' field. In this case we do find a 
systematic dependence on box size (results change by of order 20 % when we 
double the box size) so the results should be treated with some caution. For 
this case we normalised to COBE and computed the degree scale anisotropy 
using exactly the same physical resolution. 

In this paper we extend the techniques developed in [5] for calculating the 
large angular scale anisotropy, including the effects due to photons scattering 
off electrons in the ionised medium. The effect of photon drag on the velocity 
of the ionised electrons (and baryons) is a strong function of redshift, the ratio 
of the collisional damping rate f^^ to the expansion rate being given by 
t/td ~ [(1 + Z)/206j ^h^^. Since the redshift of last scattering is somewhat 
smaller than 200 in the theories we are interested in, we make the approximation 
of ignoring photon drag, so that the velocity of the scattering electrons simply 
equals the local velocity of the dark matter. This simplifies the calculation 
enormously - we do not need to evolve the baryon-electron-photon fluid. A 
useful check on this approximation is that for scale invariant adiabatic initial 
conditions we should reproduce the results for an (unphysical) standard CDM 
model without recombination. We have found reasonable agreement for the 
spectrum of multipoles obtained using our technique and that of Crittenden et. 
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al. [21] obtained using the Boltzmann equation [22]. 

Furthermore, for reasonable values of the current cosmological parameters 
at the relevant epochs the universe was matter dominated and close to critica 
density. This means that effectively only one parameter enters our calculation; 
- the horizon scale when optical depth was unity. One can simply rescale thi 
angular size of our maps to obtain results for other values of the cosmologies 
parameters Cl, fls, A, and the ionisation fraction /. 

As in [5], we solve the linearised Einstein equations analytically by Fourie 
transforming them and decomposing them into scalar, vector and tensor com 
ponents. The relevant metric perturbations are obtained as analytic integral 
of the appropriate pieces of the scalar field stress tensor, computed using FFT'i 
at each timestep. We perform a Monte Carlo calculation of the temperatun 
anisotropy by simulating photon trajectories and computing the energy shif 
each photon experiences directly. One must choose a large enough number o 
photons so that the 'shot noise' due to random scatterings, particularly thi 
Doppler term from late scatterings, is averaged out on the angular scales o 
interest. In practice we have found that for a few hundred photons along eacl 
line of sight was suflacient to suppress the white noise component in the multi 
pole spectra down to angular scales correponding to one or two grid spacings 
A typical run employs a box of 100^ for the field evolution, 50^ for the metrii 
perturbations and photon trajectories, with 200 photons along each line of sigh 
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to make a 50^ map. It takes around 5 hours on the Convex C-220 at Princeton. 
Our technique may be readily extended to compute polarisation effects. 

In synchronous gauge the temperature anisotropy experienced by photons 
along any path segment is given to first order by the Sachs Wolfe formula 

6T ^ 

where n is the direction vector of the photon, hij{r],yi) is the perturbation in 
the spatial metric, and r] is conformal time {ho^ = in synchronous gauge). 
In synchronous gauge the velocity on the dark matter is zero at all times (in 
this linear regime of course), and there is no explicit Doppler term. However, 
after integration by parts the expression may be decomposed into the following 
pieces 

ST 1 f 

— (n) = - - d,7(/ig. + h^o + hjjfi) iv, nv)nV + [-n.Vxiv, x)] f (2) 

where the surface term represents the Doppler effect of the photon scattering 
off moving electrons. This last term was not included in our previous work [5] , 
as it is unimportant on angular scales larger than the horizon at last scattering. 
In Fourier space the expression for the velocity potential x is 

XM =8nGj\v'{j^m' - ^V-'v" - ^r?-V)(e + 2es)(r?',k) (3) 

where O, and Qs are scalar components of the source fields stress tensor. Ex- 
plicit expressions for the h^J^'^ are given in [5]. 
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The initial conditions in these theories are taken to be zero temperatun 
distortion and zero metric perturbation. Each of the physical quantities o 
interest in the present calculation - the path dependent gravitational redshif 
term, the Doppler term, and the 'intrinsic' anisotropy term, can be shown to bi 
dominated at the appropriate horizon crossing epoch. Thus the calculation i 
insensitive to the early behaviour of the source fields, and issues of 'compensa 
tion' [23] do not arise. To construct photon trajectories we evolve each photoi 
backwards in time along paths starting at the position of the observer. We com 
pute the probability of a scattering event occurring during each timestep, an( 
then decide whether one actually occurred by calling a random number. Thi 
information needed to reconstruct the entire photon trajectory (initial positioi 
and direction, plus an appropriate random number seed) are stored and late 
used to step along the photon paths during a simulation. 

For the earliest part of its trajectory, a photon executes a random wall 
with step length much less than the horizon scale. During this early phase 
under the assumption that the metric perturbation is constant over the scale; 
traversed by the photon, it can be seen from (1) that the photon aquires ai 
energy shift of just —^h. We use this as an initial condition at the start o 
each photon trajectory, at the time when the photon mean free path equal; 
one grid unit. We compute the scalar 'velocity potential' x aiid the metrii 
perturbations hfj^''^ in (2) in real space as the fields evolve, and use them t( 
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determine the energy shift experienced along each photon path segment, and 
at each scattering event. 

Results 

Maps of the CMB anisotropy produced by global strings {N = 2), 
monopoles (iV = 3), textures {N = 4) and a representative case of 'nontopolog- 
ical texture' {N = 6) are shown in Figure 1 a-d. Figure le shows an equivalent 
map for the standard CDM theory [24]. The maps have been smoothed with 
a Gaussian window of FWHM 1.5°, and the angular scale subtended by the 
box side is G = 30° and 20° for the field ordering theories and standard CDM 
repectively. These angular scales are for h = 0.5, fis = 0.0125, / = 1, O = 1, 
A = 0. For other values of the cosmological parameters fi, h, f^B, A, and the 
ionisation fraction / the maps can be simply rescaled, for example if f2 = 1, 
A = the angular sizes scale as ^ a (Os//i5o/0.05)3. 

In Figure 2, we show the multipole spectra for A'' = 2, 3, 4 and 6, computed 
by averaging over the six maps obtained by observers looking at each face of a 
cube. At small I the results are consistent with previous all sky simulations [5] . 
There is significant scatter ('cosmic variance'), particularly at lower I, which 
probably accounts for all the apparent 'features' on the curves - none appear 
statistically significant. There is a striking suppression of power at higher I, 
most significant in the texture and nontopological texture theories. At I ~ 100, 
there is ten times less power in the texture theory than in standard CDM, if both 
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are normalised to COBE. Note however that this holds only if = 1, A = 0. Ii 
an open universe, or in a A dominated flat universe, there would be suppressioi 
of the anisotropy on COBE scales, requiring a larger value for e, and increasinj 
the power on degree scales. 

Is the suppression entirely a result of reionisation? We have comparec 
these spectra with those produced by standard CDM with complete reionisatioi 
assumed, calculated both using our techniques and also from the Boltzmam 
equation [21]. For that case, /^c; is approximately flat, but with a small Dopple 
bump, up to around I ~ 70, before falling off to very small values for I > 100 
The lack of high I power in the texture and nontopological texture cases seem; 
to be a result both of reionisation, and the fact that the anisotropy come; 
mainly from a time dependent potential with a large coherence length. 

In spite of the strong suppression of higher I moments, there is some powe 
up to very large I. Some of this high I power is undoubtedly attributable to thi 
localised defects themselves, but the best way to look for defects is certainb 
not in the power spectrum. 
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NonGaussianity 

The simplest measures of nonGaussianity are the third and fourth irre- 
ducible moments of the temperature distribution, known in dimensionless form 
as the skewness and kurtosis. Using the six maps produced for each theory, 
we have calculated the mean skewness and kurtosis. None of these appeared 
to be statistically significant - i.e. the mean never deviated by many standard 
deviations from zero. 

Another crude test is the level of maxima and minima in the maps com- 
pared to the r.m.s. value. The monopole maps show strong positive and nega- 
tive peaks at a level of ±3.5c7. The texture maps also appear to be significantly 
nonGaussian - the cold spot in Figure Ic is a 3.7cr event, and other maps showed 
similar hot and cold spots. It is worth emphasising that strong cold spots with 
peak level around {6T/T) ~ 5 x 10~^ are expected in the texture theory, and 
are hard to mimic with astrophysical sources. The N = 6 maps appear more 
Gaussian, in line with expectation, and have maxima and minima at around 
the 2.5cr level, similar to the standard CDM map. These numbers are sugges- 
tive, but a more detailed assessment of the level of nonGaussianity requires 
many more realisations. We have applied an optimised statistic for detecting 
nonGaussianity [25] to scans simulating the Schuster et. al. experiment on our 
maps, and do find the nonGaussianity to be detectable with signal to noise ra- 
tios as low as 1.5. It is clear however that detecting the nonGaussianity requires 
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a sophisticated statistical treatment, as well as very clean experimental data. 

Conclusions 

We have developed the techniques required for a full relativistic computa 
tion of the degree scale CMB anisotropy in nonGaussian theories of structun 
formation based on symmetry breaking and phase ordering. It is now straight 
forward to construct a statistical ensemble of microwave sky maps from whicl 
results for any experiment may be predicted, and confidence levels deduced 
Our main conclusion is that the goal of distinguishing theories based on fiek 
ordering from Gaussian theories based on inflation should be readily achiev 
able from microwave anisotropy data available in the near future. One strikinj 
difference is in the spectra of multipoles, caused partly by the reionisation ex 
pected in theories with topological defects. But even more distinctive is thi 
nonGaussianity of the anisotropy pattern - which may be detectable by thi 
degree scale experiments currently in progress. 
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Table 1. Standard Deviation on 1.5° scale 
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Table 1. The standard deviation of the anisotropy maps after smoothing wit 
= 1.5°. All theories have been normalised to fit the standard deviation or 
by COBE, {5T/T)ioo = 1.1 ± 0.18 x IQ-^. The errors quoted on the fiel 
estimates of our (one sigma equivalent) systematic errors. Statistical errors 
number for standard CDM is the value for the map shown in Figure le [1^ 
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Figure Captions 

Temperature anisotropy maps in the global string, monopolc, texture, 
and standard CDM scenarios. The side of the box subtends an angt 
field ordering cases, 20° for standard CDM. 5T /T is given in units of t 
ordering theories for N > 2, where 0o is the symmetry brciaking scale 
where /i is the string tension. Normalising to COBE. c = 5.9 ± 1.1 x 10 
1.6 ± 0.3 x 10"'' for iV = 3^ 4 and 6 respectively while SttG/^ = 5.8 ± 2. 
All the maps have been smoothed with a Gaussian of FWHM 1.5°, i 
the South Pole experiments of refs. [13], [15]. 

Power spectra for the degree scale temperature anisotropy in the g] 
texture, nontopological texture and standard CDM scenarios. If the ' 
is expanded as a sum of spherical harmonics, {5T/T) = T,aijnYi„-,{6, Ci 
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